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ABSTRACT 

Many microgravity space-science experiments require vibratory acceleration levels unachievable without active isolation. The Boeing 
Corporation’s Active Rack Isolation System (ARIS) employs a novel combination of magnetic actuation and mechanical linkages, to address these 
isolation requirements on the International Space Station (ISS). ARIS provides isolation at the rack (International Standard Payload Rack, or ISPR) 
level. 

Effective model-based vibration isolation requires (1) an appropriate isolation device, (2) an adequate dynamic (i.e., mathematical) model of 
that isolator, and (3) a suitable, corresponding controller. ARIS provides the ISS response to the first requirement. This paper presents one response 
to the second, in a state-space framework intended to facilitate an optimal-controls approach to the third. The authors use "Kane’s Dynamics” to 
develop an state-space, analytical (algebraic) set of linearized equations of motion for ARIS. 
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ABSTRACT 

Many microgravity space-science experiments require vibratory 
acceleration levels unachievable without active isolation. The Boeing 
Corporation’s Active Rack Isolation System (ARIS) employs a novel 
combination of magnetic actuation and mechanical linkages, to 
address these isolation requirements on the International Space Station 
(ISS). ARIS provides isolation at the rack (International Standard 
Payload Rack, or ISPR) level. 

Effective model-based vibration isolation requires (1) an 
appropriate isolation device, (2) an adequate dynamic (i.e., 
mathematical) model of that isolator, and (3) a suitable, corresponding 
controller. ARIS provides the ISS response to the first requirement. 
This paper presents one response to the second, in a state-space 
framework intended to facilitate an optimal-controls approach to the 
third. The authors use “Kane’s Dynamics” to develop an state-space, 
analytical (algebraic) set of linearized equations of motion for ARIS. 

INTRODUCTION 

The vibratory acceleration levels currently achievable, without 
isolation on manned space structures, exceed those required by many 
space-science experiments (DeLombard, et al., 1997; NASA 
Specification Number SSP41000, Rev. D., 1995; Del Basso, 1996; 
Nelson, 1991). Various active isolation devices have been built to 
address this need. The first in space was called STABLE 
("Suppression of Transient Accelerations By LEvitation"), which uses 
six independently-controlled Lorentz actuators to levitate and isolate at 
the experiment (or sub-experiment) level (Edberg, et al., 1996). It was 
successfully flight-tested on STS-73 (USML-02) in October 1995. 
Marshall Space Flight Center (MSFC) is developing a second- 
generation experiment-level isolation system (g-LIMIT: “GLovebox 
Integrated Microgravity Isolation Technology”), building on the 
technology developed for STABLE (Whorton, 1998). This 
compact system will isolate microgravity payloads in the 
Microgravity Science Glovebox (MSG). 
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A second experiment-level isolation system, called MIM 
("Microgravity Vibration Isolation Mount"), was launched in the 
Priroda laboratory module which docked with Mir in April 1996 
(Hampton, et al., 1997). MIM uses eight Lorentz actuators, with 
centralized control. It has supported several materials science 
experiments since its implementation in May 1996. A modified 
version of MIM (MIM II) supported additional experiments on STS- 
85 in August 1997. 

Boeing’s Active Rack Isolation System (ARIS), in contrast to the 
above payload-isolation systems, has been designed to isolate at the 
rack level; an entire International Standard Payload Rack (ISPR) will 
be isolated by each copy of ARIS on ISS. The Risk-Mitigation 
Experiment (RME) for ARIS was conducted in September 1996 on 
STS-79 (Bushnell, 1996). Each of ARIS’ eight electromechanical 
actuators requires a two rigid-body model; when the ISPR (“fiotor”) is 
included, the total isolation-system model contains 17 rigid bodies. 

In order to provide effective model-based isolation, the task of 
controller design requires prior development of an adequate dynamic 
(i.e., mathematical) model of the isolation system. This paper presents 
a dynamic model of ARIS, in a state-space framework intended to 
facilitate design of an optimal controller. The chosen approach is the 
method of Thomas R. Kane (“Kane’s method”) (Kane and Levinson, 
1995); the result is a state-space, analytical (algebraic) set of linearized 
equations of motion for ARIS. 

THE CHOICE OF KANE’S METHOD 

There are fundamentally two avenues for deriving system 
dynamical equations of motion: vector methods and energy methods. 
Both avenues lead to scalar equations, but they have different starting 
points. Vector methods begin with vector equations proceeding from 
Newton’s Laws of Motion; and energy methods, with scalar energy 
expressions. The former category uses approaches built around (1) 
Momentum Principles, (2) D’Alembert’s Principle, or (3) Kane’s 
Method; and the latter, (1) Hamilton’s Canonical Equations, (2) the 
Boltzmann-Hamel Equations, (3) the Gibbs Equations, or (4) 
Lagrange’s Equations. 
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Although some problems might lend themselves better to solution 
by other approaches, Kane’s method appears in general to be distinctly 
advantageous for complex problems. As a rule, of the above 
approaches, those that lead to the simplest and most intuitive 
dynamical equations are the Gibbs Equations and Kane’s Equations. 
And of those two approaches the latter is the more systematic and 
requires less labor. The reduction of labor is particularly evident when 
one seeks linearized equations of motion, as proved to be necessary in 
the present case (due to the otherwise excessive algebraic burden). 

An overview of Kane's approach to developing linearized 
equations of motion is presented in (Hampton, et al., 1998), along with 
a summary of the relative advantages of the method. See Kane and 
Levinson (1979, 1985) for more extended treatments. 

DESCRIPTION OF ARIS 

The total dynamical system S consists of the stator S (ISS and 
the integral frame, from the motion of which ARIS isolates the ISPR), 
the flotor F (the ISPR), eight electromechanical actuator assemblies, 
and the umbilicals. (See Figure 1.) The flotor is connected to the 
stator by the eight actuator assemblies, and by a variable number of 
umbilicals. The actuator assemblies also (and fundamentally) act as 
the vibration isolation devices. 

Each actuator assembly consists of a Lorentz (voice-coil) 
actuator, an arm, an upper stinger, a push-rod, a lower stinger, and a 
position sensor. (See Fig. 2 for a kinematic diagram, and Fig. 3 for a 
CAD drawing, of a single actuator.) One end of each actuator arm is 
connected to the flotor through a cross-flexure which allows the flotor 
a single rotational degree of freedom with respect to the stator. The 
other end of the arm is connected to one end of the push-rod through 
the upper stinger, a wire of very high torsional stiffness. Each upper 
stinger provides two rotational degrees of freedom in bending. The 
opposite end of the push-rod is connected to the stator through the 
lower stinger, another short wire which allows three rotational degrees 
of freedom (two in bending, one in torsion) with respect to S. 

Each stinger is modeled as a massless spring. The umbilicals are also 
considered to be massless; they are modeled together as a single, 
parallel spring-and-damper arrangement, attached at opposite ends to 
stator and flotor at effective umbilical attachment points S u and F u , 
respectively. This effective umbilical applies both a force and a 
moment to the flotor. The force is assumed to act at point F u . 



Figure 1. ARIS Control Assembly 


The stator, the flotor, and each actuator arm and push-rod are 
considered to be rigid bodies, with mass centers at points 
S\ F\ A*, and P*. respectively. The superscript * indicates the mass 

center of the indicated rigid body; the subscript i corresponds to the i th 
actuator. (/ = 1,...,8). All springs (cross-flexures and stingers) are 

assumed to be relaxed when the ISPR is centered in its rattlespace (the 
“home position”). 

COORDINATE SYSTEMS 

With the ISPR in the home position, fix eight right-handed, 
orthogonal coordinate systems in the flotor, one at each of the cross- 
flexure centers. Let the i ,h coordinate system (/ = 1,...,8) have origin 

Fj (/ = 1,...,8) located at the center of the i th cross-flexure, with axis 
directions determined by an orthonormal set of unit vectors 
f [j = 1,2,3). (The overhat indicates unit length, the index i 
corresponds to the i h actuator assembly, and the index j distinguishes 
the three vectors.) Orient the unit vectors such that f ^ is along 

the i ,h arm, toward the i th voice coil; / is directed parallel to the other 

segment of the i Ih arm and toward the upper stinger (which is located 
~ i *i *i 

at Aj ); and / is in the direction f ^ X f ^ (along the intersection of 
the two cross-pieces of the i th cross-flexure). 


l.ORUNTZ COIL. 



Figure 2. Kinematic Diagram, Including the i’ h Actuator 
Assembly and the Umbilical 

Fix a similar right-handed coordinate system 1*2,3) in 

the arm of each actuator. Locate each system a - such that it is 

coincident with the corresponding flotor-fixed coordinate system 
~ i 

f when the flotor is in the home position. 

At the respective lower stingers (points S,), place eight push-rod- 
fixed coordinate systems p l , , and eight stator-fixed coordinate 
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Arm 


Figure 3. A Single ARIS Actuator 

systems £*• . Orient these 24 coordinate systems such that when the 
stingers are relaxed (i.e., with the ISPR in the home position), the 
coordinate directions p\ , and £^ are co-aligned for the i !h actuator, 

with p ' (along with s 2 , in the home position) directed from 5/ 
toward A 2 . 

Define finally a primary, central, flotor-fixed, reference coordinate 
system with coordinate directions f . All other flotor-fixed coordinate 

systems / . are assumed capable of being referenced (e.g., by known 
direction cosine angles) to this system. [See Equation (4).] 

ROTATION MATRICES 

Let the a_ } coordinate system rotate, relative to the / coordinate 

^ i 

system, through positive angle q\ about the / axis. Similarly, let the 
orientation of the aj coordinate system, relative to the p, coordinate 

system, be described by consecutive positive rotations q l 2 (about the 
axis) and q 3 (about the moved 3-axis). And let the orientation of 

the p . coordinate system, relative to the iV coordinate system, be 

described by consecutive positive rotations q l 4 (about the £3 axis), 
q l 5 (about the moved 2-axis), and q' 6 (about the moved 1-axis). 

Let c l j and s'j represent the cosines and sines of the respective 

angles q l j . Then the rotation matrices among the several coordinate 
systems for the i th actuator assembly are as follows: 

e]\ 44 

■ P 2 ‘= - 44 +444 44 
»' 44 + 444 ~ C 4 S 

(_3 J L 



("If -Af 

S 4 C 5 

-•% £1 


i i 

C 5 S 6 1^2 

6+444 

/ i 

c 5 c 6j[£ 3 



and 




( 2 ) 


(3) 


Finally, define a rotation matrix between the eight flotor-fixed 

a i 

coordinate systems / and the single, flotor-fixed, reference 


coordinate system / 




fn 

fn 

/2I 

f 22 

/23 

fi \ 

fn 

/*33 



(4) 


GENERALIZED COORDINATES FOR S 

The 48 angles q' } are the generalized coordinates of the system. 

For the / th actuator the six associated generalized coordinates are as 
follows: q[ is the angle at the cross-flexure of the i lh actuator; and 

q‘ } are the angles at the upper stinger; and q\ , q\ and q‘ 6 are the 
angles at the lower stinger. 

GENERALIZED SPEEDS FOR S 

Define generalized speeds U } for the system as the time rate of 

change of the generalized coordinates of S in the inertial reference 
frame: My = q l j (for j— /=1 8) . (5) 

ANGULAR VELOCITIES OF REFERENCE 
FRAMES AND RIGID BODIES 

Designate the reference frames corresponding to the stator, the 
/' lh push rod, the i lh arm, and the flotor, by the symbols S , P { , A, , 

and F , respectively. Let 5- and F i represent, respectively, the 
coordinate systems in S and F defined respectively by 

[fj £2 £3] and / f ^ f ^ . Two intermediate reference 

frames were introduced previously to permit describing the angular 
velocity of each push rod relative to the stator; designate those 

intermediate frames corresponding to the / th actuator assembly by 
Rj and Q { . Another intermediate reference frame was previously 

introduced between frames P t and A, ; designate this by T i . 

Let each intermediate reference frame have a frame-fixed, dextral 
set of unit vectors. Indicate the unit vectors for each of these frame- 
fixed coordinate systems by using the corresponding lower case letter 

( r*j corresponding to /?, , etc.). The following, then, give the 
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expressions for the angular velocities of the various reference frames 
and rigid bodies of S : 


F ‘^=u\r v Ti <* Ai =“& 

(6- 8) 

S ‘ (ti R ‘ =M4£ v R ‘ CD^ =1*5 £2, and Q ‘ CO*' = u‘ 6 q\ 

( 9 - 11 ) 

Using the addition theorem for angular velocities, the angular 

velocities of the rigid bodies of S are 


5f 0) A ' =U l 2 p 1 +M4I3 +« 5£2 

(12) 

5, CO / ‘ = M4J3 +« 5£2 + t 

(13) 

and V ” u \p\ + u ?Lt> + u \Li + ~ li \ f ^ 

(14) 


LINEARIZED ACCELERATIONS OF THE CENTERS OF 
MASS FOR THE RIGID BODIES OF S 


Taking the time derivatives of the respective linearized velocity 
vectors yields expressions for the linearized accelerations of the 
centers of mass, for each rigid body. Note that the linearized velocity 
vectors may be used in this step — the full nonlinear accelerations need 
not be determined. This is a tremendous savings of effort, which 
would not be afforded if Newton's Second Law were applied directly, 
instead of Kane's approach. 

(26) 




S )a A ' = [— 44 - (a‘ 2 + + (a,'4 + 


c 


2 U 2 


A \ u 5 


(a' 2 +l' 2 )u 6 }p r 


(27) 


BASIC ASSUMPTIONS 

In the subsequent development of the ARIS equations of motion, 
it is assumed that ARIS works as intended; i.e., that the ARIS 
controller prevents the ISPR from exceeding its rattlespace constraints. 
It is also assumed that the small-angle approximations hold for angles 

q l j . Angular velocities and angular accelerations are assumed to be 

small as well. This means that the use of first-order linear 
perturbations will permit the full nonlinear equations of motion to be 
approximated accurately by a set of first-order linear differential 
equations. Finally, it is assumed that the angular velocity of the stator 
is negligible, and that the stator translational velocities and 
accelerations are small. 


LINEARIZED VELOCITIES OF THE CENTERS OF MASS 
FOR THE RIGID BODIES OF S 

Represent by r AB the position vector from arbitrary point A to 


arbitrary point B. Define the following position vectors, using the 
indicated scalars: 




: A f £2 +/ 2&. 


r t * 


■4K 


S,P 

r" =p 2 p 


—2 


(15-17) 


r** = a;a; + a\i S'j , and /' F = 'f + f { /, . (18. 1 9) 


First time derivatives of the appropriate position vectors, under 
the stated assumptions, yield expressions for the velocities of the 
centers of mass, for the seventeen rigid bodies. The following 
expressions are the linearized velocities for those centers of mass. 
(The pre-subscript indicates that the expressions are linearized; the 
pre-superscript indicates the reference frame assumed fixed for 
purposes of the differentiations.) 


iV 

- 

r4*i +4iJ 

(<= L. 

...8); 



(20) 

V* 


M-k +/ik 

Jf|+I 

r . . 

44 + 

44)4, 




*u 

■) U 2 ip 2 

+ 4 )‘t 




(21) 

and s \i 

{ 4 

fW\ ~44 -( V 2 +4)'L +V.x]p] 




♦i 


- V '('4 +«(,) + V l(“' 

+«i)J P 


-v|«4+(v2+4 

klfy 

(22) 

where 

v{ = 

V* 2 

=/i- 

-4 

and = /q 

(23 

-25) 


and s \a = [/j'u , 1 -v^ -(v 2 +l$u\ + 44 ]// 

+ [-/iV “4 (4 + 4>)+ v '(4 + 

+ [vj«2 - + ( v 2 + 4 )«6 }py < 28 ) 


LINEARIZED PARTIAL VELOCITY VECTORS FOR THE 
POINTS OF S AT WHICH THE CONTACT/DISTANCE 
FORCES ARE ASSUMED TO ACT 

The partial velocities and partial angular velocities are formed by 
inspection of the relevant velocity vectors. These partial velocities are 
then (and the order here is crucial) linearized by neglecting higher 
order terms. 


LINEARIZED PARTIAL VELOCITIES OF P* 

For the r lh push-rod the linearized partial velocities are 


S \£' =P (for r = 1,2,3), =-p'‘j'£\+q\i^\ 


(29, 30) 


S, P, i i M 

iZs =P 2^6lr 


and 


Si Pi i ( 

l —6 =P 2 V 


-till +i;)(3i. 32 ) 


LINEARIZED PARTIAL VELOCITIES OF A? 

For the / lh arm the linearized partial velocities are 


=Q, =- a WiP 2 + (°2 + a Wi)F i ' 

(33, 34) 

= -(a ‘ 2 +a\q\)p y +(b,' “44 + a\q' 2 F y 

(35) 

=44 +44 + 4 )/> +(4 -44)^2 _ (44 +44 +4 i)p 2 • 

(36) 

= (44 +44 +l We)p ] -44p' 2 +(44 -4 )p y » 

(37) 

s ‘, Vs = -44 £ 2 + (4 + 44 + l i )p y 

(38) 


LINEARIZED PARTIAL VELOCITIES OF F* 


For the flotor, the linearized partial velocities are 


S , F 


-1 


-1 


* -/i7 v 


(39) 
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V'2 =-( V 2'/2 +V \)t[ + ( V 1<?3 + v 2 ~ V Wl ~fW\ )Py (40) 

=[-^<h -V 2 +/!4)p\ +(vl -v^ +A4)p 2 +4<&Pf < 4 1 > 

s ',u =[- v k' - v 2 + v 'k +4}~f'4 -4kj 

+ (v] - + vl,2/5 + f\q\ - (v|<7^ + v\q 5 + 1 3275 j/J^ , (42) 

vr =kk +4) +4 + i 'dlp\ -4 &p 2 -k +4d +A4)p\ - (« > 

and V{ =k V 2^ +V ' K + 4d + V 2 ~4fk ~f\4 +4 p\- (44) 

LINEARIZED PARTIAL VELOCITIES OF F„ 


F F 

Define measure numbers for r “ as follows: 

f r --Xr£ + r r.£, + z '£ 


Since 


5, F u s. F . 
\Vr=lYr + 


3 f s -(0 F xr F F - 

du r 


(45) 

(46) 


the linearized partial velocities for the umbilical attachment point 
F u can be expressed as follows: 


' 5| vf”=Vf +W,-W 2 ’ '>2"=>2 - Z f»P 2 +Y ^' (47 ’ 48 ) 

5 i^=V<‘ -y Fu ^F U p 2 . K“= 5 ' ( _wf-z fu p; + x fu £;. <49, 50) 


-l 


5, F u S, F 


S { . F„ _5, . F 
l 


+Z F uP-X F «P, • S \£J)£ ~Z Fu P 2+ Y Fu P 3 . (51,52) 


LINEARIZED PARTIAL VELOCITIES OF £ 

The linearized partial velocities of F l are 

S i vf- =0, s }v£ =-fe +493)23. ^=^-4^. (53-55) 

■V4 = fe -4k -fe +4k+k +4k +%)+4^kr (56) 

•V? = [- /; U + 4 )+ 4^6 k; + 4/2 + ^ k' 2 + (4 - 4<4 k > (57) 


and S ‘l—6 = 4<?2 £2 - (4 + l Wi - 4 k 3 • 


(58) 


LINEARIZED PARTIAL ANGULAR VELOCITIES 
FOR THE RIGID BODIES OF S 

The following are the linearized partial angular velocities for the 
system. 

For the / th pushrod: =0 (for r = 1, 2, 3) , (59) 


s ‘i =iv 


S, ,P t m 
I —5 “ —2 ’ 


and ^5)5 = <7*. (60-62) 


For the « th actuator arm: =0, ^O)^' = p‘ , 


(63, 64) 


5 ja£' =[3, =|3, ^6)5' = £2 , and ^CO^' - gj (65-68) 


5. A A 

' fD 1 — c 1 fi'k ' — y 


S: A: 


For the flotor: 5 j tof 


iv 5 /“2 =p|’ 


-1 s, 


-1 


CO 3 =£.3,(69-71) 


1^4 —^3 ’ ;®5 - Ll ’ 1—6 *£ 


(72-74) 


and — 0 ( r — 1 6; I — 2,..., 8) , (75) 

LINEARIZED ANGULAR ACCELERATIONS 
FOR THE RIGID BODIES OF S 

LINEARIZED ANGULAR ACCELERATION 
OF ACTUATOR PUSH-ROD P, 

S )a P ’ =u‘ 6 p i +u‘<,P 2 + “4/4 • (76) 

LINEARIZED ANGULAR ACCELERATION 
OF ACTUATOR ARM A ■ 

' S ' ( ot A ' = («2 + Me,kl + «5£*2 + (rij + i<ik 3 • (77) 

LINEARIZED ANGULAR ACCELERATION OF THE FLOTOR F 



CONTRIBUTIONS TO THE SET OF GENERALIZED ACTIVE 
FORCES DUE TO THE RIGID BODIES OF S 
CONTRIBUTIONS DUE TO THE FLOTOR 

On orbit (i.e., neglecting the effects of gravity), the flotor is acted 
upon by forces and moments due to the each Lorentz coil, actuator 
arm, and umbilical; and by direct disturbances, 
c c 

Let — F ' and — M_ ' represent, respectively, the force and 
moment exerted by the / lh Lorentz coil (located at A \ ) on the flotor, 
where jF C ’ = F C ' a J , assumed to act at point F ( , (79) 

and M_ C ‘ = r F,A ’ x£ C ' = -F C ‘ (/( + /' )a‘ . (80) 

Let F 1 and M_ ' represent, respectively, the force and moment 
exerted by the I th actuator arm on the flotor, at the I th cross- flexure. 
Since F^' is a noncontributing force, it can be ignored in the analysis. 

p 

The total moment M_ due to the eight cross-flexure springs has 

value M F =tk‘q\ \f (81) 

I _ ' 

where k[ is the / ,th cross-flexure spring stiffness. 

Let F U and M_ U represent, respectively, the force and moment 
applied to the flotor by the umbilical; where the force is assumed to act 
at flotor-fixed point F u . Umbilical force F_ U is given by the equation 

F U = (- V, - c l-*1 )4, + (- k 2 x 2 - C 2 X 2 )l 2 

+ (- £ jc 3 - c 3 x 3 )| 3 + F _ h ; (82) 

where ^ is some appropriate stator-fixed coordinate system; x x , x 2 , 

and x 3 are the umbilical elongations in the respective (i = L2,3) 
directions; F_ h is the umbilical bias force in the home position; 
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k 2 , and £ 3 are umbilical spring stiffnesses; and c x , c 2 , and c 3 are 
umbilical damping constants. Umbilical moment M_ U is given by 

a/ /= (->M>i -Yi<h)| +(-ic 2 < t>2 -1^% +M h , (83) 

where § { , (j) 9 , and <|) 3 are components of the umbilical angle of twist 
<|> in the respective ^ (i = 1 ,2,3) directions; A/^ is the umbilical bias 

moment in the home position; K ( , K 2 , and K 3 are torsional 
umbilical spring stiffnesses; and Yj , y 2 , and y 3 are torsional 
umbilical damping constants. 

Let F°and M_ D represent, respectively, the unknown 
disturbance force and moment acting on the flotor. Assume F_ D to act 
through the flotor mass center F*. Define F t D and M ^ to be the i th 

components, respectively, of F_ D and M_ D , componentiated in Fy . 

In terms of the above, the flotor’ s contribution to the set of 
generalized active forces, for the I th generalized speed, is 


,#=V ?A v J f+yc\ q \ 1 m 


F ? 

The umbilical contributions to the ; Q r ’s , viz., 

V" A s ;o/ / , are addressed in the following two sections 
The remaining terms of the l Q , ’s are as follows. 

-/!>?■ 

V +4 

+«+*,v -«• 

v +/?(<« +A-fu\ 

Vi .(-^1=0, 

ivf fM-vy -^+/,y)^+(vi -m +/MH’. 


(85, 86) 
(87, 88) 

(89) 

(90) 
(91,92) 

(93, 94) 


s \o£ • f i /^/ 3 -*iV + ^iV +*iV 

‘ l — \ 


(95) 
(96, 97) 

(98) 

(99) 


V £“=[-'! -vW V +*fo +4)-/,* 

^ y)+v!-vy ■ 1 ifK ' ;i 


+ f3® 


X • m d =U -*< +< , (loo 

X s 'd\-bf i] h s )d (-^)- o° 2 - 103 > 


lU 

X- 


r s 


i^/ 3 |=x- x*m/ 3 . 


(104) 


/I’f •£ D =^ +V 2U + 4) +V 3] / )° 

+N(<fc+?6)+ v 'b -4^2 +N . 005) 

X ■ M r> =(</,' + A/f -(4 +^j,< , (106) 

;vT(-^ : )=^[/;4+fc-4kL 00?) 


=0, (108,109) 


)^^=^(i[ + iU-4 x(^?7; 

:, ,vf -f 0 =V^' -$? +fe +7 +V& +v!, ,(110) 

;a('M"= l ,a( /W D , X ■(-f')=-f C 'l i d, (111, 112) 

]o( (-A/ ; )=0, and X VX' 


013, 114) 

V J v 

Notice the coupling between the unknown-disturbance measure 
numbers and the generalized coordinates. This coupling will make the 
disturbance input matrix E [in Eq. (203)] time-varying. 

UMBILICAL FORCE F U 

Equation (82) expresses umbilical force F^ in terms of 
umbilical-elongation components Xy , x 2 , and jr 3 , and their time 
derivatives. These items must be re-expressed in terms of the 
generalized coordinates and generalized speeds. 

If the umbilical attachment point F u is at F uh in the home 
position, then 

Xj = ( r SuF “ -r s ‘ F “ h )• | , for / = 1, 2, 3 . (115) 

But r S “ F “ -r S ‘ F ‘* =r F ^ =r S ' F ' +r F ' F ‘ -r V “ -r S ^ , (11 6a, b) 

where the right-hand-side terms can be expressed by 

L' F ' =l\p 2 -l\a x -l\a 2 , r F ' F ‘ =x Fu f+y F J 2 +Z Fu f 3 (117,118) 

r 1 “ ~Vu^.i Xvul "l" Xv;l $ ■ an 0 L “ “ h = X 1 _ l"3'o£2 ^"^ 0^3 • C 1 9, 1 20) 
for appropriately defined coefficients. 

Define now the following rotation matrix: 
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r \] 

r 12 r 13 _ 


-i] 

£i 

I, 

> = 

r 2\ 

r 22 r 23 


^1 

£2 

-1 



/31 

r 32 r 33_ 


£3 


( 121 ) 


In terms of the coordinate system, Eq. (116b) can now be written ; 


r FuhF " =*i4, +x 2 t, 2 +x& 3 . 


where 


for 




^1 

r \2 ^3 


C, +y F L 1 -‘h~‘i*\ +z FA ~ 1 Wa +/ ‘U + ^) 

“ / \ /it 1 ^ 1/1 i\ 

*1 


r 2l 

r 22 r 23 


c> -Z K W -q\ + 44 ) 

3j 


3\ 

r 32 6.1 


C 3 + (4 - ks + “ ^1 + ^6 ) + (V /6 


Ci - x F - x s 


(122) 

(123) 

(124) 


C 2 = y Fu - y Su + 1\ -/, - y 0 , and C 3 = z Fu ~Z Su ~ Zo- (125,1 26) 
Differentiating, 
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V 


1 

r-j 

<-*> 

1 


*2 

> = 

r 2\ r 12 r 23 


A 


Jm hi h3_ 



>'F U (“! -«3 ~ l 3 U 4 +/ l(«3 +M l)j 

~ Z F ( w 2 +H 6.) - '' r F u ( M I ~ u \ +u 4 


[4 ~ X F U M +(»„ “4 1'4 +u f) + h u 6 


(127) 


Using Eqs. (82) and (123) through (127), a linearized expression could 
now be written straightforwardly for the umbilical force F_ U . 

UMBILICAL MOMENT 

Equation (83) expresses umbilical moment M_ L in terms ot 
angle-of-twist components <|>j , <J) 2 , and <j) 3 , and their time 

derivatives. These items must be re-expressed in terms of the 
generalized coordinates and generalized speeds. 

Let represent the rotation of the flotor, relative to the stator, 

from the home position, n ^ is the rotation axis, and <|> is the angle of 
twist about that axis. Note that <|) f = (|)rtx • for / = 1,2,3 . ( 1 28) 


-1 -1 -1 

Express as + £ 2 /, + g 3 /, • 


Define rotation matrix Q by 


r -1 ' 

/, 


f - 1 1 

s . 

- 1 
V-2 

•=lef 

-i 

£2 

{l\ 


-1 

£3 


(129) 


(130) 


The linearized 3x3 rotation matrix ( Q has elements / Qj; defined as 


follows: 


f -1 1 

L 

r 

- 1 

\L 

> = 

[ l \ 

- 


l l l 

<7l — <73 —^4 

1 

Q 5 


-q\ +<?3 + ^i 


-^2-96 


1 


+ 

1 


<72 +<76 


-1 


-1 


(131) 


(132) 


For small <|> it can be shown (Salcudean, 1991) that 
0 -g 3 g 2 

<l> g3 0 ~g\ 

_~gi g\ 0 _ 
where the post-superscript T indicates matrix transposition and tr , Q 
represents the trace of { Q . Substitution from Eq. (131) into Eq. 


iQ-i <l_ 

(l+tr ,qT ’ 


(129), and simplification, yields g t = 

0 

\q\ -(?*)■ 


and g 3 


1 


033) 
(134, 135) 


1 l 

82 = ~~"T <7s ’ 

<P V 

Substituting from Eqs. (133)-(135) into Eq. (129), and transforming 

into the coordinate system by use of ( Q , one obtains the following 

expression for the spin axis: 


/0 


Since has unit length, 


Ad-Ji-q if 


\n 


(136) 


(137) 


Use of Eqs. (128), and (136) leads to the following linearized forms 
for angular position and rotation rate: 


and 


<t>l 


'll 

r l2 

<3 




► = 

*21 

hi 

h3 


<75 

03 J 


/31 

hi 

h3_ 


</i ~<?4 

y 


'll 

’ll 

<3 


~[ u i +*4) 

02 

► = 

h\ 

hi 

r 23 



03 


h\ 

hi 

<33. 


Mj M 3 — M^ 


(138) 


039) 


From Eqs. (84), (138), and (139), a linearized expression could 

now be written straightforwardly for the umbilical moment M_ U . The 

flotor’ s contribution to the set of generalized active forces, for the r* h 
generalized speed, could then be found by substituting the expressions 

for F [ U (previous section) and M U , into Eq. (84). 


CONTRIBUTIONS DUE TO THE ACTUATOR ARMS 

The forces and moments acting on the t lh actuator arm are due to 
the respective Lorentz coil (located at A 3 ), the flotor (through the 
/ lh cross-flexure), and the respective push-rod (through the upper 
stinger). The coil force F_ ' is the only contributing force. The 
contributing loads, in the above indicated order, are as follows: 


F C ‘ = F C> a j , assumed to act at point F- , 

(140) 

M q = r F ' K x£ c - = (/,' + l‘ 4 )a 2 x F C ‘ q\ = -F C ’ (lj + 1 

(141) 

-M_ F ' = —k[q\ h , and =-k 2 q 2 p } 

(142, 143) 

where and l\ are pertinent geometric lengths, and k 

2 and /: 3 are 

pertinent upper-stinger spring stiffnesses. 


In terms of the above, the contribution for the / actuator arm to 

the set of generalized active forces, for the generalized speed, is 

,Qr' -= S i)lA ■ K C '+ S I \k a ' +K c ' -K f '} 

1. (144) 

The individual terms of the ( Q ’s are as follows: 


■ F_ C ‘ =0, S [(of (m a + M_ C ‘ - M_ F> )= 0 , (145, 146) 

s t v Fi ■ F_ Ci =0, S ,c^ [m* +^' -M Fi )=-kirf 2 ’ 

(147, 148) 

S iY.3 Z C ’ =F C 'l[ , 

(149) 

^ (m A +K C ‘ -M F ' )=k{qi —k\q' 3 — F C ' (/' 

+ /i), (150) 

s lV j-F c '^F c '(n-iil 

(151) 

5 ,cof +m Ci (/; +/;) 

1. (152) 

S ,v F '-F C '=-F c Ali( q ' 2+q ‘ 6 )-l^}^ 

(153) 

S ,u A (m a +M Cl -M F, )=-F C, [ l \ +l < 154 ) 

]v *' 2 -F C ' =0, and (a / 4 -M Fi )=~k 'A - 055, 156) 
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Notice the coupling between the control inputs and the generalized 
coordinates. This coupling will make the disturbance input matrix B 
[in Eq. (204)] time-varying. 

CONTRIBUTIONS DUE TO THE PUSH-RODS 

p 

The contributing loads on each push-rod are moments M_ ' and 
— M_ A ‘ , where (using pertinent lower-stinger stiffnesses) 

M 1 ’ =<4a 1, -k'd? 2 -Kit, 1 ■ 057 ) 

The contribution for the /' lh push-rod to the set of generalized 


CONTRIBUTIONS TO THE SET OF GENERALIZED INERTIA 
FORCES DUE TO THE RIGID BODIES OF S 

Represent by / A k 1 A ' the central moment/product of inertia of the 

i'* actuator arm for the j and k body-fixed coordinate directions 

^ * p j p* 

a- and a* k . Define push-rod inertias If ' analogously, where the 
single subscript indicates that the axes are assumed to be principal 

F i F* 

axes. Let I - k represent the central inertia scalar of the flotor for 

the flotor-fixed coordinate directions f . and f ^ . Use the symbol H 

to represent an angular momentum vector. Associated post- 
superscripts on H have the same meanings as for the inertias. The 

contributions to the generalized inertia forces for S can now be 
expressed. 

CONTRIBUTIONS DUE TO THE PUSH-RODS 

The contributions (, Q * ) ' to the generalized inertia forces due to 


the / lh push-rod are as follows: 
(for r - 1,2, 3;/ = 


pjp- 


(163) 


'-,h pM ' 


(idf = s, i if (-ifA y-i? lPi its . 


(164) 

(165) 


Ur = 5 ;vf V : } 


S ,.P, 
6 


~,K 


pjp; 


( 166 ) 


,Qr‘ XX \M P ' 

P , 

The individual terms of the { Q/ s are as follows: 

(158) 

Then tl 


= 0 (r = 1, 2, 3) , 

(159) 

1Y.2 

X s ■fe''’-"*) 

^3^3 — ^4^4 • 

(160) 

s i ( X- 

X^-aC') 

■— in 
■— m 

1 

11 

(161) 

S, A* 
^3 

and X ' -M. A ) 

'=k' 2 q' 2 -k‘ 6 q‘ 6 . 

(162) 


CONTRIBUTIONS DUE TO THE ACTUATOR ARMS 

The contribution (/ 0* ) to the generalized inertia forces due to 
the i lh actuator arm, for the r lh generalized speed, is 

j. (167) 

{[Q r ) ’ *S are as follows: 


5, A 
- m A, i“ 


1^2 ^ -ly}^ (r/3 (169) 

uthur 


UF+ 




(170) 


(171) 


s, a' ( S, / 
I —4 I ~ m A, ,9 


(«i) 2+ a 2{ a 2 +l 7i) 

Uj+(4 + iif 


s w A ' 

l —4 


S, A, 


-iK 


A, I A- 




{ _3 /L L 


,/C M ' 


(172) 

(173) 




S ,.A 


Vc’ • -m, ',a 


S, A, 


A, /“ 


v r 

■/ . \ 

? . n 

r" w n 

(«{) 

S=. 

Ut ~ 
1 


— a \ a 2 U 2 — ( fl 2 "6 ^3 ) “6 } > 


(174) 


• -,K 


aA/a' 


= -4k fe +«4)-4k <4 -4k’ («2 +«4). (i75) 


/L 6 


S, 4 ( 5 A. . 

'tv,' -I -m A \a 1 | =~m A 


l Aj li 


“2 

( ^2 + ^3 F “6 “ a l ( + ^3 )^5 

s , <o? ■( -tii 


and 5 ,0>M- I- 


(176) 

(177) 


CONTRIBUTIONS DUE TO THE FLOTOR 

The contribution (/(?*) to the generalized inertia forces, for 
the r lh generalized speed, is 


(iQ'X -\£ {- m F 5 )9 YdUs: 

V J V 


,fIf- 


(178) 


8 



Then the individual terms of the (/£?*) s are 




S /Z2 j = -»i F {v]/]'« 

Rv!)MvO : 


|«2 — V } V 3 li 3 -vjVjui 


-V|V2«5 + 


(viF +vj(vi+/j) 






i f* r & f* 
& I - m F l“ 


= -Wr 


-/M-/M 


t/| l -vjv^2 + 


(179) 

(180) 


(181) 

(182) 


4 +( V 'Fl«3+|"( V lF +V 2( V 2+4) 


^ -V2V3M5 -vjv_X 


]<£{-, if 


FI F* *i ,F//* 


'32 V 


33 


(183) 
fe+^-lijj, (184) 


iiU j--m / . |-/2 (/{ V* 2 )h,' 


EQUATIONS OF MOTION FOR THE SYSTEM 
KINEMATICAL EQUATIONS 

There are 48 kinematical equations for the system, one for each 
generalized speed: Uj = q } (for j = / = 1,...,8) . (191) 


DYNAMICAL EQUATIONS 


Six dynamical equations are obtained using the following 
process. First, add the respective contributions of the 17 rigid bodies 
to the set of holonomic generalized active and holonomic generalized 
inertia forces, for each generalized speed (i.e., r=I,...,48). The 
holonomic generalized active force for the 4 generalized speed is 

F r = l{F r )> , (192) 

> = 1 


where ( F r ) J is the contribution to the set of holonomic generalized 
active forces due to the f h rigid body. That is, 

F=,g+ htf+htf. ( 193 ) 

i=l /=! 

Likewise the contribution to the set of holonomic generalized inertia 
forces is 



) j ={,q:Y + i{ l Q;f 4 {,q;T 


(194) 


Second, develop the relationship between the dependent and the 
independent generalized speeds in the form: 

u) = i,A rs ii\+B r (i = 2 8; j = r = 7,.. .,48) . (195) 

. 1=1 

Where in the above equation u] are the six independent generalized 



(185) 

(186) 


(187) 



(188) 


(189) 


(190) 


speeds, and Wy are the 42 dependent generalized speeds. A rs is a 

42x6 matrix, derived from the nonholonomic constraint equations 
(see next section). The nonholonomic and holonomic generalized 
active forces are related to each other as follows: 

~ 48 

F r = F r + ^ F s A sr (r = 1,...,6) . (196) 

s =7 

Similarly, the nonholonomic and holonomic generalized inertial forces 
are related to each other as follows: 


48 


f; =?; + if; A sr (r = 1 6). 

. t =7 

(197) 

Kane’s Dynamical Equations, then, are F r +F r = 0 , 

(198) 

for r = 1,...,6 . 


CONSTRAINT EQUATIONS 



The kinematical and dynamical equations together are Fifty-four 
in number: 48 kinematical, 6 dynamical. Since the complete set of 
equations for an n-degree-of- freedom system numbers 2n, and since 

the system S has 48 degrees of freedom, 42 more equations are 
needed to describe completely the motion of the system. These 
missing equations are the holonomic constraint equations (in 
nonholonomic form) for the dependent generalized speeds 

u ) (i = 2,.. .,8; j = \ 6). 

Since the velocity and angular velocity of the flotor center of 
mass F* is the same irrespective of the actuator path chosen for 
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describing its position, a set of constraint equations can be written in 
vector form using the following: 


s, r 


fs -dr s ' r ' 


dt 


(=*;/' 0 = 1, 7 = 2... .,8), (199) 


and 


s \o/' = \(0 F ‘ (i = 1,7 = 2 8). 


( 200 ) 


If one expands Eqs. (199) and resolves them into a common 
coordinate system (here, the / coordinate system), one obtains the 
following twenty-one (motion) constraint equations: 
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2,.. 

.,8; 

j = 1.2.3). 

(201) 


Similarly, if one expands Eqs. (200) and resolves them into a common 
coordinate system (here again, the /. coordinate system), one obtains 
the remaining twenty-one (motion) constraint equations: 
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ki 

fh 



fh 
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«3 

fh 



fh 


“5 

fh 


kJ 


-fh 
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fh 


u 2 

fh 
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fh 


u 4 

fh 


«5 

fh 




,8; j = 1,2,3). (202) 


(It should be noted here that, although Eqs. (201, 202) are in 
nonholonomic form, the constraints they represent are actually 
geometric.) 


STATE-SPACE FORM OF THE EQUATIONS OF MOTION 

The dynamical, kinematical, and constraint equations can be 
arranged in the following state-space form: 


"/ o' 

f q | 

'O N 

1-1 

'o' 

t i 

' O' 


\ ” = 




\ij+ 


O M 

l“ , J 

K C 

vr 

B 


E 


where the state vector consists of the 48 coordinates q and the 6 

independent generalized speeds u f , associated with actuator #1; the 

constant submatrices M, K , and C are system mass, stiffness, and 
damping matrices, respectively; the symbols l and O represent, 
respectively, an identity matrix and a zero matrix of appropriate 
dimensions; vector i contains the control currents to the Lorentz coils; 

and vector d_ is the disturbance vector. The input matrices B and E 
are time-varying matrix functions of the coordinates. N is a constant 


submatrix that incorporates the kinematical equations and the 
holonomic constraints. 

The disturbance term [e] {d} accounts for the umbilical bias 
force F^ h and moment M _ h , and the unknown direct disturbance force 

F° and moment . Recall that in the development of the 

foregoing equations the angular acceleration of the stator was assumed 
to be negligible. However, the translational acceleration of the stator, 
although presumably unknown, ca mot be neglected. In fact, that 
acceleration is the source of the umbilical contribution to tlotor g- 
jitter. To include this indirect disturbance contribution, one simply 

adds an (unknown) indirect acceleration disturbance term a 1 to each 
of the Eqs. (26), (27), and (28). Along with the other disturbances, 
this indirect disturbance will appear in the final term of Eq. (203). 

MODEL VALIDATION 

AUTOLEV software, marketed by Online Dynamics, Inc., was 
used to create a full nonlinear model of ARIS, including the actuator 
(rigid-body) dynamics. AUTOLEV was then used to develop and 
verify the linearized equations presented in this paper. 

The linearized model (without umbilicals) was implemented in 
MATLAB m-files, using actual values of system parameters. The 
mathematical model was then checked for kinematical consistency. 
The procedure used was first to compare the eight position vectors 
from a common point on the stator to the flotor center of mass, as 
traced through the eight actuators, with the flotor centered in its home 
position. The eight position vectors matched exactly. The procedure 
was repeated with the flotor moved from its home position, in six 
degrees of freedom. The position vectors tracked within acceptable 
limits. The linearized MATLAB model was compared with the full 
nonlinear AUTOLEV model. The static response of the linearized 
model tracked the static response of the nonlinear model, for small 
angles. 

An independent model was developed using the DENEB 
Envision software, with current CAD models of an ARIS-outfitted 
ISPR. This model was used as an independent (static) check of the 
actuator kinematics. 

FUTURE WORK 

The next tasks will be the addition of umbilical forces to the 
AUTOLEV and MATLAB equations, and the dynamical validation of 
the MATLAB model. The procedure will involve simulating the 
application of various loads to the flotor, and verifying that the eight 
position vectors track in this dynamic simulation. System dynamics 
will be incorporated into the Envision model, along with the capability 
of state-space discrete-time control. The MATLAB and Envision 
models will then be available, respectively, for centralized, state- 
space/optimal controller design and for closed-loop system simulation. 
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